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1 .  INTRODUCTION 

Recently  R.  F.  Link  (1981)  used  finite  difference  operators  for  a 
unified  approach  to  determining  the  moments  of  discrete  probability  dis¬ 
tributions  . 

He  presented  derivations  for  several  well-known  distributions:  the 
binomial,  Poisson,  geometric,  and  hypergeometric.  Following  his  approach 
we  give  derivations  for  several  important  distributions  in  occupancy  prob¬ 
lems  and  some  other  well-known  discrete  distributions  (negative  binomial, 
power  series,  and  factorial  series) .  Also,  Johnson  and  Kotz  (1981)  gave 
the  formulae  of  (descending  or  ascending)  factorial  moments  for  the  binomial, 
Poisson,  hypergeometric,  negative  binomial,  and  negative  hypergeometric. 

If  the  moments  about  the  origin  are  expressed  by  means  of  finite  dif¬ 
ference  operators,  the  corresponding  descending  factorial  moments  can  be 
found  immediately  by  a  theorem  given  below. 

In  this  article,  E  is  the  displacement  operator  and  A  is  the  dif¬ 
ference  operator;  i.e. 

Ef(x)  =  f(x+l)  and  Af (x)  =  f(x+l)  -f(x). 

* 

This  work  was  done  while  the  author  was  visiting  the  Department  of  Epidemi¬ 
ology  and  Public  Health  at  Yale  University  in  1981  and  was  supported  by 
Academia  Sinica. 


The  relationship  between  E  and  A  is: 

A  =  (E— 1)  ;  E  = (1+A) . 
If  X  is  a  discrete  random  variable  and 


Pj.  =  P(X=j)  j  -  0,  1,  2,  ... 


th 


Its  s  moment  about  the  origin,  y'  is 

V  ■  Vo  -j  J- 


and  its  sth  descending  factorial  moments  is 

fs") 

where  j  =  j  (j-1)  ...  (j-s+1)  . 

In  what  follows, =  0,  for  a<b  or  b<0  ,  and  An0m  = 


,n  mi 
A  x 


x=0 


2.  THE  CLASSICAL  OCCUPANCY  PROBLEMS 


The  basic  models  of  the  classical  occupancy  problems  are  as  follows: 
there  are  m  urns  and  n  balls,  and  the  balls  are  distributed  among  the  m 
urns.  Let  be  the  number  of  urns  containing  exactly  t  balls  (t=0,  1, 
...,  n) .  Find  the  distributions  of  M  (cf.  Johnson  and  Kotz,  1977). 

According  to  whether  the  urns  and  balls  are  distinguishable  (DT)  or 
indistinguishable  (IDT),  and  whether  empty  urns  are  permitted  or  not,  there 
are  several  possible  models,  four  of  which  are  listed  in  Table  1,  where 
T{i}  (n’m)  denotes  the  number  of  ways  in  which  n  balls  are  distributed 
among  m  urns  under  the  model  {i}  (cf.  Johnson  and  Kotz,  1977,  p.37). 
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Table  1.  The  Classical  Occupancy  Models 


Model 

Urns 

Balls 

Are  empty  urns  permitted 

T{i}  (n,  m) 

{1} 

DT 

IDT 

Yes 

/  n  +  m  -1  j 

\  m  -1/ 

{2} 

DT 

IDT 

No 

U:i) 

{3} 

DT 

DT 

Yes 

n 

m 

{4} 

DT 

DT 

No 

Am  0n 

The  models  {1}  and  {2}  correspond  to  the  Bose  -  Einstein  system  and 
the  models  {3}  and  {4}  correspond  to  Maxwell  -  Boltzmann  system.  Let  M  ^ 

be  Mt  under  the  model  { i }  (i=l,  2,  3,  4).  It  is  well-known  that  the  distrib¬ 
utions  of  M  are: 


(-l)j 


m-r  w  n-Hn-  (  r+j  )  ( t+1)  -1 
j  A  m-r-j-1 


{2} 

P{Mt  =  r}  = 


C-Dj 


j  m-r\,  n-(r+j)t-l\ 
l  j  A  m-r-j-1  1 


PlM^. 


{3} 


r}  = 


(") 


n 


m 


m-r 

£ 

j=0 


(-1)' 


m-r 


n : 


5  /  (t!)r+j(n-(r+j)t)! 
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{4} 


(“) 


P{M  =  r}  = 

t  ,m  „n 


m-r  .  /  m-r  \ 

!  <-DJ  (  3  ) 


n ; 


A  0  j=0 

We  can  express  these  in  a  unified  form: 


(t!)r+j(n-(r+j)t)! 


Am-r-j  Qn-(r+j)t 


{i}  (“) 

P{M  =  r}  =  r 


L{ 


1  m-t 

1  r(n>^y  z  ("DJ  T{i}  (n-(r+j ) t ,m-r-j )  h±  (n,t,r+j) 
l}  *  j=°  for  1=1,2, 3, 4, 


where 


h^(n,t,a)  =  1*2  (n,t,a)  =  1 


h  (n,t,a)  =  h  (n, t , a)  =  - - - : -  • 

J  4  (t!)a(n-at) ! 

Let  be  the  sth  moment  of  11  about  the  orgin.  We  can  derive  by 

ills  t  &  ills  7 

using  finite  difference  operators.  From  the  definition  of  moments,  we  have 

m  {i}  m  r - i 

p,:,  =  E  rS  p{M.  =  r}  =  E  P{M^^  =  r}  Er  0S 

{l}s  r-0  C  r-1 

{i} 

and  we  can  express  P{Mt  =  r}  as 


■f*T  ^  i  *  j1,  i  i  m 

T{i}  (n,m)  =  r}  =  E  (-l)r  J  1/  \  T{±}  (n-jt,m-j)  h.(n,t,j)  . 

j=r 


Thus  we  get 


:{i} 

(n,m)  u 

{il 

s 

1 

m  m 

E  E 

:=1  j=r 

(-1) 

r+j 

(ffl: 

[u>  (n-^- 

>m-j  ) 

hi(n,t,j)  Er 

m 

>  E 

j=l 

(-l)j  | 

iml 
u , 

ITU} 

(n-jt. 

m-j) 

h. 

l 

(n,t,j) 

i  (-Dr 

r=l 

it ) 

Er  0S 

m 

I  m\ 

‘  E 

j=l 

(-l)j  | 

L! 

ITU} 

(n-j  t , 

m-j) 

h. 

l 

(n,t,j) 

(1-E)^  0S 

s 

■  E 

j=l 

U  Tm  <n-Jt’ 

m-j)  h 

±  (n. 

t,j' 

1  Aj  0S, 

for 

i=l.  2,  3,  4, 

4 


i. e.  , 


m 


(2-l)  ^  |jj  T{±j  (n-jt,  m-j)  h±  (n,t,j)  AJ  0S/T{i}  (n,m) 


3.  THE  RESTRICTED  OCCUPANCY  PROBLEMS 

In  the  models  above,  if  the  maximum  number  of  balls  allowed  in  any  urn.  is 
k,  we  get  the  restricted  occupancy  problems.  The  models  are  as  follows: 

There  are  m  urns  and  n  balls.  Each  urn  contains  k  cells.  The  n  balls  are 
assigned  among  the  m  urns  so  that  each  cell  contains  at  most  one  ball. 

Fang  (1982)  considers  several  models  according  to  whether  urns,  balls  and 
cells  are  DT  or  IDT,  and  whether  empty  urns  are  permitted  or  not.  Ten  such  models 
are  listed  in  Table  2. 

Fang  (1982)  obtained  the  distributions  of  Mt  under  the  model  (i),  i=l,  2, 

.  ..,  10.  They  have  the  following  unified  form  except  for  model  (9): 


P{Mt=r} 


T^(n,m,k) 


m-r 

E 

j=0 


(n-(r+j ) t  , 


m-(r+j),  k)  gi(n,k,t,r+j) 


where  T^(n,m,k)  is  the  number  of  ways  of  distributing  n  balls  among  m  urns 
under  the  model  (i)  and 


1 

for 

i  =  1,  2, 

(*>a 

for 

co 

II 

•H 

(k)0  n  (at) 

for 

i  =  5,  6, 

n! 

for 

i  =  7,  8, 

(t !) a(n-at) ! 

n! 

for 

i  =  10. 

m(a)(t:)a(n-at) ! 
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Table  2.  Models  of  the  Restricted  Occupancy  Problems 


Symbol  if  empty 
urns  are  permitted 

Symbol  if  empty 
urns  are  not  permitted 

urns 

balls 

cells 

(1) 

(2) 

DT 

IDT 

IDT 

(3) 

(4) 

DT 

IDT 

DT 

(5) 

(6) 

DT 

DT 

DT 

(7) 

(8) 

DT  ' 

DT 

IDT 

(9) 

(10) 

IDT 

DT 

IDT 

We  see  that  the  distributions  of  have  a  form  which  is  similar  to  that 
of  the  classical  models.  Thus 


s  rm  i 

O.D  y-(^  -  E  [jJ  T.  (n-jt,  m-j,k)  g.(n,k,t,  j)  aV/T.  (n,m,k)  i»l,2 . 8.10  , 

s  J  t 


where  =  E  (Mfc)  under  the  model  (i) 


4 .  COMMITTEE  PROBLEMS 

A  group  contains  n  individuals,  any  w^^  of  whom  can  be  selected  at 
random  to  form  the  ith  committee  (i=l,  2,  ...,  r) .  Find  the  probability 
that  exactly  m  individuals  will  be  committee  members.  This  has  been 
called  the  committee  problem  which  has  been  studies  by  several  authors 
(cf.  Johnson  and  Kotz,(1977),  and  Holst  (1980) ) . 

Let  L  be  the  number  of  individuals  who  do  not  serve  on  any  of  the 
r  committees;  then  exactly  (n-L)  individuals  will  be  committee  members. 
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We  can  find  the  distribution  of  L  by  using  the  finite  difference  operator 
(White,  1971)  or  by  the  inclusion-exclusion  principle  (Sprott,  1969).  It 


is 


n-k 

P{L=k}  =  E  (-1)3 
j  =  0 

x  th 

Its  s  moment  y  is 


r 

n 

i=l 


n-max  w .  n-k 
(4.1)  y;  =  E  1  E 
k=l  j=0 


(-dj  ncji  w//» 


i=l  '  w . 

x 


w. 

i' 


J  v 

-  ”  ”  \)  Ek0S 

i-i  v  »i  /  wy/ 


Ek  0S 


k=l  ct=k 


In  particular,  if  w.= 


w  =  w,  we  have 
r 


(4.2) 


5.  OTHER  DISTRIBUTIONS 

(1)  Negative  Binomial  Distribution 

The  formula  for  the  moments  of  the  negative  binomial  distribution 
as  obtained  by  Link  (1981)  contains  double  summations,  and  Link  suggests 
that  a  way  should  be  found  to  get  rid  of  these  . 
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We  can  express  the  negative  binomial  distribution  as 


(Patel,  Kapadia  and  Owen,  1976,  p.22).  Now 

<5-«  s-'j-ora  *kY°s 

=  pk  (1-qE) ~k  0S 
=  pk  (p-qA)  _k  0S 
=  (1  -  q/pA)  k  0S 

(2)  Power  Series  Distribution 

Here  p.  =  a  0j  /  f (0)  j=0,  1,  2,  ... 

J  J 

Where  f ( 0)  =  E.°°  a  0^ 

J=0  j  . 

The  class  of  power  series  distributions  includes  the  Poisson,  binomial, 
negative  binomial  distributions.  From  the  definition,  we  have 


(5.2) 


y" 

s 


=  x  a.0^E^OS/  f (9) 

j=l  J 


=  2  a .  0J  2  |J\A1 

j=l 


i=0  \  i 


5  E  2  a.0J*  (  ~j  \  A* 

i=l  j  =i  d  *  1  / 


'  ?  -£  fC1)  i1  0s 
i=l  x‘ 


os  /  f (0)  . 

9s  /  fee) 

/  f(0) 
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(0)  denotes  the  ith  derivative. 
Factorial  Series  Distributions 


where  f 
(3) 


(i) 


0(j)  Ajf(0) 

Pj  j !  f (6)  j  -  0.  1>  2 s  . . . 


e(J)  j 


where  f  (0)  2  — r1  A  f(0)  and  f(0)  satisfies  suitable  conditi 

j=0 

(cf.  Johnson  and  Kotz,  1977,  pp. 87-88).  We  have 

(5.3)  f ( 0)  u'  =  E  -4-H —  A^  f  (0)  0S 

j=0 


ions 


,(j) 


=  E 
j=0 


rr -  AJ  f(0)  E  (J  )  A1  p 


i=0 


o  CO  -f  o 

=  E  E  0(j)  AJ  f (0)  -fS- 
i=l  j=i  ( j -i) !  ~  1* 

_  ft  r  e(a+i)  .oH-i  rrm  A1  0s 

-  z  £  — Zl -  A  f  (0)  -  »■  -- 

•  1  «  *  1* 
1=1  a=0 


=  E  0(l)  A1  E  ^  Aa  f  (p)  AJL 

•  -i  0i  •  1  • 

1=1  a=0 


s  /  •  \  •  .1  „s 

=  e  0  1  a1  f(e-i) 
i=l  1* 


6.  DESCENDING  FACTORIAL  MOMENTS 

The  formulas  which  have  been  found  by  Link,  Johnson  and  Kotz,  and 
by  us  can  all  be  written  in  the  following  form: 
u'  is 

s  =  E  a.  A  0  .  s  =  1,  2,  ... 

i=l  1 
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It  is  well  known  (cf.  Johnson  and  Kotz,  1977,  pp.  8-10)  that: 


(6.1)  Xs  =  E .  1  X(J)  Aj  0S  /  j! 

J=1  J  » 


(s) 


thus  there  is  a  relationship  between  and  y  '  which  is 

s 

(6.2)  \i'  -  E  p(j)  Aj  0S  /  j! 

S  i=l 

Let  A  =  (a„)  be  a  s  x  s  matrix,  where 


u.  -f  AV/j! 
l  0 


for  i  >  j 
for  i  <  j . 


Obviously,  A  is  a  lower  triangular  matrix  with  a^  =  1  (i=l,2, .  . .  ,s)  .  Hence 

A  is  non-degenerate.  From  this  fact,  we  obtain  the  following  theorem. 

Theorem  If  the  moments  about  the  origin  have  the  form 
s 

y'  =  I  a.  AJ  0S  s  =  i  o 

s  -iJ  ~  s 

J=1 

then  the  corresponding  descending  factorial  moments  are 

y^  =  a  s !  s=l,  2,... 

s 

By  using  the  theorem  and  formulas  obtained  by  Link,  Johnson  and  Kotz,  and 
us,  we  immediately  derive  the  descending  factorial  moments  of  the  follow¬ 
ing  distributions.  These  moments  can  be  found  in  textbooks  or  papers,  but 
the  derivations  are  often  elaborate  or  lengthy. 

(1)  Binomial  Distribution 

(s)  _  (s)  s 

y  -  n  p  . 

(2)  Poisson  Distribution 

(s)  s 

y  =  t 


(3)  Geometric  Distrbution 


,(s) 


=  s.  q 


s-1 


/  p‘ 
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(4)  Hypergeometric  Distribution 
U(S)-„<S)  X<S)  ./  N(s). 

(5)  Negative  Binomial  Distribution 

[s] 


.Cb)  = 


(s) 


-  (k+s-1)  (q/p)  =  k  (q/p) 


[s] 


where  k  =  k  (k+1)  . . .  (k+s-1) 


(6)  Power  Series  Distributions 
,(s)  _  as  .(s) 


y 


=  0  f 


(0)  /  f(6) 


(7)  Factorial  Series  Distributions 

/S>  =  0(s)  AS  f  (0-s)  /  f(0)  . 


(8)  The  Classical  Occupancy  Problems 
,(s)  _  (s) 


y  =  m 
{i} 


T{i}  (n  st,  ms)  h^n.t.s)  /  (n,m)  ,  i=l,2,3,4, 


(9)  The  Restricted  Occupancy  Problems 
C  s)  (s) 

y  =m  T.  (n-st,  m-s,  k)  g.(n,k,t,s)  /  T.  (n,m,k)  i=l,2, . . . ,8,10. 

(i)  1  1  x 


Fang  (1981)  obtained  the  same  results  by  using  characteristic  functions. 


CIO)  Committee  Problems 
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